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Abstract 

Entanglement is defined for each vector subspace of tlie tensor 
product of two finite-dimensional Hilbert spaces, by applying the no- 
tion of operator entanglement to the projection operator onto that 
subspace. The operator Schmidt decomposition of the projection op- 
erator defines a string of Schmidt coefficients for each subspace, and 
this string is assumed to characterize its entanglement, so that a first 
subspace is more entangled than a second, if the Schmidt string of 
the second majorizes the Schmidt string of the first. The idea is ap- 
plied to the antisymmetric and symmetric tensor products of a finite- 
dimensional Hilbert space with itself, and also to the tensor product 
of an angular momentum j with a spin 1/2. When adapted to the 
subspaces of states of the nonrelativistic hydrogen atom with definite 
total angular momentum (orbital plus spin), within the space of bound 
states with a given total energy, this leads to a complete ordering of 
those subspaces by their Schmidt strings. 



1 Introduction 

If quantum entanglement jl| is to be regarded as a physical resource ^ , then 
it seems sensible to consider the entanglement not only of individual states, 
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but also of collections of states of a given composite quantum system. On 
the other hand, it is not clear how to combine measures of entanglement of 
individual states in such a collection, because of the possibility of superpos- 
ing given states to form new ones. Since the natural organizational unit for 
any collection of states is the vector (sub) space spanned by those states, we 
are led to the problem of quantifying the degree of entanglement inherent 
in a given vector subspace of the whole state space of a quantum system. 
Examples of vector subspaces of interest might be the space of states with a 
given total energy, or the space with a given total angular momentum. More 
generally, a state subspace might be labelled by the eigenvalues of any incom- 
plete set of commuting observables. A situation that arises often is one where 
the state space is associated with a tensor product representation of some 
symmetry group or algebra, and the subspaces of interest carry irreducible 
subrepresentations of that algebra or group. Again, the example of angular 
momentum springs to mind; we may be interested in a subspace of states 
carrying a definite total angular momentum, for a quantum system made up 
of several subsystems, each contributing angular momentum to the total. In 
such cases, the problem of quantifying the entanglement of individual irre- 
ducible subspaces is seen to have an essentially group-theoretical character; 
any measure of entanglement of such a subspace must surely involve such 
group-theoretical constructs as the Clebsch-Gordan coefficients of the corre- 
sponding group or algebra. Conversely, considerations of the entanglement 
of such irreducible subspaces seems likely to throw interesting new light on 
familiar group-theoretical reduction problems. 

In what follows, we consider only bipartite systems, and vector subspaces 
V of a complex, finite-dimensional state space H = TCi ®7i2, where the two 
factor spaces have dimensions di and d2 respectively, and are equipped with 
the usual scalar products. Extensions to the multipartite case seem hkely to 
face the same sort of difficulties as entanglement measures for state vectors 
of multipartite systems. 

Example 1 : As more specific motivation, consider the following 3-dimensional 
vector subspaces of H, in the case where Hi = H2 has orthonormal basis 
{ei, 62, 63, . . . e^v}, > 3: 

(A) Va is spanned by the orthonormal vectors (ci 62 — 62 ® ei)/\/2, 
(62 (8) 63 - 63 62) /V^, (63 (g) ei - ei (8) e3)/V2. 

(S) Vs is spanned by the orthonormal vectors (ei 62 + 62 ® ei)/\/2, 
(ei (8) ei - 62 <8) 62) /V^, (ei (8) ei 62 (8) 62) /V^. 

The bases of Va and Vs so defined consist of three maximally entangled 
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vectors in each case. However, it is possible to find another basis in Vs with 
only one entangled vector, namely the set {ei®e2+e2(8>ei)/v^, CiCSci, 62062}, 
whereas every choice of basis in Va consists entirely of maximally entangled 
vectors. It is intuitive that Va is 'more entangled' than V^, and we seek to 
quantify such differences. 

2 Operator entanglement and subspaces 

In general, any vector subspace V <H = T-Ci®H2, of dimension I < d < did2, 
is characterized by a corresponding hermitian projection operator P, 

PV = V, P^ = P = P\ Tr(P)=d, (1) 

and we suggest that measures of entanglement of the operator P provide 
suitable measures of entanglement of V. 

Measures of operator entanglement have been considered previously in 
other contexts |31 E]. The central idea is to consider each linear operator 
(matrix) A as an 'operator vector' \A)) in the ((ii(i2)^-dimensional 'operator 
vector space' E-^ of all linear operators on Ti, with Hilbert-Schmidt scalar 
product 

{{A,B)) = TriA^B). (2) 

Similarly, for r = 1, 2, linear operators on Tir can be considered as operator 
vectors in the d^-dimensional operator vector space . li A = B <^ C, 
then \A)) = \B)) ® \C)), and A is unentangled. Otherwise, A is entangled. 
The operator vector 

1^)) = ^ 1^)) (3) 

corresponding to the projector P divided by Vd, is a unit operator vector in 
E-j-c, according to ((H) and 0- We can define measures of entanglement of this 
unit operator vector in E-^, just as we define measures of entanglement of 
unit vectors in Ti. To this end, we note firstly that \P)) will have a Schmidt 
decomposition, 

|P)) = ^ \Ei)) ® + ^ IE2)) ® IF2)) + . . . \Ek)) ® \Fk)) , (4) 
where 

K < K = mm{dl, dj} , 
Pi > P2 > ■ ■ ■ > Pk > , P1+P2 + . . .+Pk = ^, (5) 
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while \Ei)), \E2)), ... \Ek)) are orthonormal operator vectors in E-j^^^, and 
I-F2)), . . . \Ek)) are orthonormal operator vectors in E-j^^. If we introduce 
the superoperator jH] density matrix, of dimension (^1(^2)^ x (^1^2)^, 

R=\P)){{P\, (6) 

and then define the reduced superoperator density matrices M*-^-* and R'-^^ 
by tracing over the second (respectively, the first) vector subspace of E-^, 
then \Ei)), I-E2)), . . • \Ek)) are eigen operator vectors of M^^-*, and \Fi)), I-F2)), 
. . . \Ex)) are eigen operator vectors of M'^^\ in each case with eigenvalues 
Pi, P2, ■■■Pk- Furthermore, the unentangled unit operator vector closest 
to \P)) - in the sense of the norm defined by the scalar product - is 
I -El)) l-^i)), and its distance from \P)) is 

^dUp))) = {m-m®{{F,\){\p))-m®\F,)))f' 

= y/2(l-Vpr). (7) 

This distance provides a partial measure of the entanglement of \P)), and of 
V, and we shall also write it as Soiy). 

The entanglement of V is fully characterized by its corresponding K- 
dimensional 'Schmidt string' 

'5(V) = (pi, P2, - --Pk, Pk+1 = 0, pk+2 = 0, . . . , p^f = 0) . (8) 

Various partial measures of entanglement can be defined in terms of the 
Schmidt string, including Sd as above. Thus the 'information' measure of 
entanglement of |P)), and hence of V, is 

Sji\P))) = Sj{V) 

= -Tr(iWlog2(i(^)) : 

K 

= -^Pa'^Og^iPa) , 

while the 'trace' measure of entanglement is 

St{\P))) = St{V) 

= l-Tr(RW2) 



= -Tr(R(2) log2(i(2)) 

(9) 



= 1 -Tr(i(2)2) 

(10) 
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A better indicator of entanglement is provided with the help of the notion 
of majorization [HI 13 IE] ■ Thus we may say that V < H is more entangled 
than W < n, with Schmidt string S{W) = (gi, ga, • • •)> if Pi < ?i AND 
Pi + P2 < Qi + <i2 AND . . ., that is to say, if SiV) is majorized by 5(W), 
which we write as SiV) -< iS(W). When SiV) -< iS(W), it can be shown 
m that Eniy) > ^d(W), ^/(V) > £/(W) and ^t(V) > SriW). But when 
neither of S{V) and iS(W) majorizes the other, some of these inequalities 
and not others may be reversed. In that situation, it is best to say only that 
V and W are differently entangled. 

The preceding two paragraphs merely paraphrase for operator (or sub- 
space) entanglement what is well-known for state entanglement. Many state- 
ments that hold true for the tensor product space of states ?i, go over to the 
tensor product space of operators Efi, without the need for new proofs. In 
what follows, we give some properties of subspace entanglement as defined, 
and then some examples of naturally arising entangled subspaces and their 
Schmidt strings. 

Property 1: If in the situation described above. Hi is embedded as a subspace 
in a larger space H[, and H2 is embedded as a subspace in a larger space Tig, 
so that 



then S{V), when V is regarded as a subspace of Ti.[ ® 0.2, differs from iS(V) 
when V is regarded as a subspace of Hi <^Ti.2, only by the addition of the 
appropriate number of zeros on the right-hand end. In this sense, our notion 
of entanglement of a subspace is stable against embeddings. 

Property 2: If V has the form Vi ® V2, where Vi < Hi and V2 < H2, 
then \P)) = I Pi)) (g) IP2)) is unentangled, and so also is V, according to our 
definition. In this case, the Schmidt string has the form 



In particular, H itself is unentangled, and S{H) has the form ()12j] . 

Property 3: If V is 1-dimensional, spanned by the unit vector \v) say, with 
Schmidt decomposition 



V <Hi®H2<H[®H'2, Hi< H'l , H2<H'2 



(11) 




(12) 



k 




(13) 
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then 

k 

P = \v){v\ = J2 K){ep\ ^ \fa){U\ . (14) 

a, (3=1 

This defines the Schmidt decomposition of \P)), whose Schmidt string then 
has as components the PaP/3, with a suitable ordering. From this it is easily 
deduced that the Schmidt string of \v) majorizes the Schmidt string of |m) 
if and only if the Schmidt string of the subspace spanned by 1^;) majorizes 
the Schmidt string of the subspace spanned by \u). This guarantees that our 
notion of entanglement of 1-dimensional subspaces is consistent with that for 
state vectors. In particular, it is also easily seen that 

k 

a, (3=1 
k 

= -2j]p„log2(pJ = 2£,(|^;)). (15) 

a=l 

Thus the information measures of entanglement of a 1-dimensional subspace 
and of any unit vector within that subspace differ only by the constant factor 
2. 

3 Antisymmetric and symmetric subspaces 

Example 2: As a generalization of (A) in Example 1 above, consider the 
'antisymmetric tensor product' space Va < Ti. = 'Hi ® H2, where Hi = H2 
has an orthonormal basis {|ei), |e2), . . . , |en)}. An orthonormal basis for Va 
is provided by the n{n — l)/2 vectors 

\eki) ^{\ek)^\ei)-\ei)^\ek))/V2, k < I, k, I e {1, 2, . . . , n}, (16) 

and the projector onto Va is then 

n 

Pa ^ Y \eki){eki\ 

k<l=l 

1 " 

= 2 X] (|efe)(efc| ® |ei)(e;| + \ei){ei\ O |efc)(efe| 

k<l=l 

- \ek){ei\ \ei){ek\ - \ei){ek\ \ek){ei\) . (17) 
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Labelling the unit operator vectors in E-^^ and Ef^^ as 

= |ei)(ei|, |2)) = |e2)(e2|, . . . , \n)) = |e„)(e„|, 
|n + 1)) = |ei)(e2|, \n + 2)) = |e2)(ei|, . . . , |3?2 - 2)) = |e„)(ei|, 
\3n - 1)) = |e2)(e3|, \3n)) = |e3)(e2|, . . . , |5?2 - 6)) = |e„)(e2|, 

1^2)) = |e„)(e„_i|, (18) 

we then have from |T7j). the unit operator vector 

where the x matrix A with matrix elements Ars takes the form 

A = B®C®C®...®C. (20) 

Here B is n x n, with all diagonal elements equal to 0, and all nondiagonal 
elements equal to 1. Each of the n{n — 1)/2 copies of C is 2 x 2, with diagonal 
elements equal to 0, and off-diagonal elements equal to —1. 

It follows from (fT^ that the matrix elements of M*^^-* in this case are just 
those of AA\ whose eigenvalues are easily calculated from (pUj) to be 



^'^-) = M;r3T)((«-i)^i-i.--ij- (21) 

where the 1 appears — 1 times. Then 1)211) is the Schmidt string for V^- It 
follows that 

£d{Va) = v/2(l- v/(^r^WM), 
Si{Va) = log2 {2n{n - 1)^/") , St{Va) = {n + l){3n - 4)/[4n{n - 1)] . (22) 

Note that, as n — > oo, £^d(Va) and £^t(Va) tend to constants, whereas 
Si{Va) ~ log2(n). 

Example 3: Consider again the space Ti as in Example 2, and let de- 
note the 'symmetric tensor product' space of dimension n{n + l)/2, with 
orthonormal basis \eki), k < I = 1, 2, . . . , n, where 

\eki) = (|efc) ® le^) + |e/) \ek))/V2, k <l, 

\ekk) = \ek) ® |efe) . (23) 
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A similar calculation to that for the antisymmetric case shows that the 
Schmidt string for is 

where again the 1 appears — 1 times. Then 

SjiVs) = log, {2n/{n + 1)^/") , St{Vs) = {n - l)(3n + 4)/[4n(n + 1)] . (25) 

The asymptotic behaviour of these quantities as — > cxd is similar to that in 
the antisymmetric case of Example 2. 

We note from ^IT\ and ()24|) that, in Example 1, the 3-dimensional an- 
tisymmetric and symmetric subspaces of the A^^-dimensional space Ti. have 
iV^-dimensional Schmidt strings 

5(V^) = ^(4, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0,...,0) 

'5(V5) = ^(9, 1, 1, 1, 0, 0,...,0), (26) 

so that iS(Va) -< <S(ys), consistent with our intuition that Va is more entan- 
gled than Vs- 

4 Coupled angular momenta 

Example 4 '■ Consider the coupling of two angular momenta, ( Ji, J,, J3) with 
spin j, and {Si, S2, S3) with spin 1/2. In this case, the full state space is 

H = nj® ni/2 = v,+i/2 © v,_i/2 , (27) 

where the space Hk, for k = j or 1/2, or Vk, for /c = j + 1/2 or j — 1/2, has 
dimension {2k + 1), and carries the corresponding irreducible representation 
of the spin Lie algebra su{2). We are interested in the entanglement of Vj±i/2, 
regarded subspace of H. 

Let |m), m = j, j — 1, . . . , —j denote the usual orthonormal basis of 
eigenstates of J3 in Hj, and let |-|-) and |— ) denote the usual orthonormal 
basis of eigenstates of 5*3 in Hi/2- Then let 

|l)) = |+)(+|, |2)) = |-)(-|, |3)) = |+)(-|, |4)) = |-)(+|, (28) 
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defining an orthonormal basis of operator vectors in -E1/2, and let 

\m,n)) = \m){n\, m, n = j, j - 1, . . . , -j , (29) 

defining an orthonormal basis of operator vectors in Ej. Then S± = Si ±182 
and 5*3, regarded as operator vectors in E1/2, take the form 

|5+)) = |3)), |5_)) = |4)), |53)) = ^(|1))-|2))). (30) 

while J± and J3, regarded as operator vectors in Ej, take the form 

I J+)) = VOWFUj - 1)) + v/(2)(2j-l)|j - 1, J - 2)) + . . . 

... + V(2j)(l)|-J + 1,-J)) 



I J_)) = v^(l)(2j)|j - 1, j)) + v/(2)(2j-l)|j - 2, J - 1)) + . . . 

... + v/(2j)(l)|-J, + 

1^3)) = jIj, j)) + (j - l)|j - 1, J - 1)) + . . . + (-j)l - J, -j)) ■ (31) 
Recall inj that the 7i operator X, defined by 

X = J+ ® 5_ + J_ O 5+ + 2 J3 ® ^3 , (32) 

takes the eigenvalue j on the subspace Vj+1/2 and the eigenvalue — (j + 1) 
on the subspace Vj-i/2- It follows that the projector from H onto Vj±i/2 is 
given by 

^'"'■^ = ±27TT(*n^'±y + i)0- (33) 

where I denotes the unit operator on H. 

Consider firstly the projector Pj+i/2- From (jHUjl . (jSH), (jH^ and (jSBI), we 
see that this operator, regarded as an operator vector on E-j-^, and normalized 
to a unit operator vector, takes the form 

+2IJ3))® 1^3)) 

+ + + |j-l,j-l)) + ... + |-j,-j)))®(|l)) + |2)))}. 

(34) 
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Here the terms on the last hne represent the operator vector corresponding 
to the operator [j + 1)/. Expression has the general form 

O'j) 4 

\Pj+l/2)) = '^A^rn,n),a\m,n)) 0\a)) , (35) 

{m,n)={-j-j) a=l 

and we wish to claculate the eigenvalues of the reduced superoperator density 
matrix 

^^'^ = E E {Arr,n),aAl^,n),p} , (36) 

(m,n)={-j-j) Q!,/3=l 

or, what is the same thing, the eigenvalues of the 4x4 matrix Q with elements 

Qafi = E {"^(m,"),a"^(m,n),/3} • (37) 

{m,n)=(-i,-i) 

The only nonzero elements are, from (j3T| and (j3l|l . 

Qn^g. - 2(/+l)(2,!l)^ ((^^- + ^)^ + (^^-)^ + - + W-) 
4j + 3 



6(2j + 1) 



(3^ 



^33 = ^44 = ^(/^ ^) \ ((l)(2j) + (2)(2j - 1) + . . . + (2j)(l) 

= ^^^^ 

and 

Qi2 = Q21 

- 2(JTT) WTI? (^^^' + ^^^^^ + ^'^^^'^ + • • • + + ^0 

2j + 3 



6(2j + 1) 



(40) 



and the eigenvalues of Q are now easily calculated to be (j + l)/(2j + 1) 
(multiplicity 1) and j7(6j + 3) (multiphcity 3). The Schmidt string of Vj+1/2 
is therefore 

+ (41) 



2j + 
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We then have as scalar partial measures of the entanglement of this subspace, 
Sd{V,^i/2) = ^2(1- v/(j + l)/(2j + l)), 

= -log2 



^|y7(2i+l) ^ l)(i+l)/(2i+l)' 



2j + l 

St{Vj+i/2) = 2j(4j + 3)/[3(2j + 1)^] . (42) 

We see that as j — >■ 00, all these quantities approach constant values. This 
is a consequence of the fact that the Schmidt string (PT|) approaches the 
constant value 

^^0 = (I, I, I, l) , (43) 



2' 6 6 6^ 

and can perhaps be understood as follows: as j gets large, the number of 
states in Vj+1/2 with larger and larger positive or negative eigenvalue of 
J3 h/2 + Ij ® S3 increases, and these states have smaller and smaller en- 
tanglement, with the entanglement reaching zero for the highest and lowest 
weight states. 

A similar calculation shows that the Schmidt string of Vj~i/2 is 

5(V,-,/2) = ^ [j, ^ j . (44) 

In this case, 

£D{V,-y2) = ^2{l-^j/{2j^), 

£:,(V,-V2) = -iog2 ( ^^^^-^ 

£T{Vj^i/2) = 2(j + l)(4j + l)/[3(2j + 1)^] . (45) 

As J — > 00, these quantities approach the same constant values as in the 
previous case. 

Note that 5(Vj+i/2) >- 5(Vj_i/2), so that Vj-1/2 is more entangled than 

5 Application: Electron spin and H-atom 



Example 5: Consider the space 'H^"'^ of bound states of the nonrelativistic 
hydrogen atom with principal quantum number n, where n is a positive 
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integer. This space is n^-dimensional, with the structure 

= e e . . . e Un-i , (46) 

where 7ii is (2/ + l)-dimensional, corresponding to the orbital angular mo- 
mentum content I — 0, 1, . . . , n — 1. Allowing for the spin of the electron, 
we have as the relevant state space including spin, 

= (Ho ® H1/2) ® (Hi ® H1/2) ® . . . ® (Un-l ® H1/2) 
= (Vl/2) © (Vl/2 ® V3/2) © (V3/2 © V5/2) © 
■ ■ • © (Vn-3/2 © Vn-1/2) ■ 

(47) 

Prom the results of Section 4, we see that the Schmidt strings corresponding 
to these subspaces are 

, 1 / , 2k-l 2k-l 2k-l 

SW = -{2k+l.— — — 

2 2 2 



~ / 2A: + 3 2A: + 3 2A: + 3 

'^^^^^ = 4(^(,'^+''^'^'^ 

f ^ 1 3 3 
for k = -,-,..., n — - . 
2' 2' ' 2 



(48) 



Now we see a remarkable ordering of these subspaces by their spin-orbit 
entanglement. From least entangled to most entangled, as indicated by their 
Schmidt strings, we have: 



12 



S(Vi/2) = (l, 0, 0, o) )- 5(V3/2) )- 



1 2ra-2 2n-2 2r2-2 



_ 1 1 1 1^ 
- '2'6'6'6'^ 



^ ~ , 1 / 2n 2n 2n\ 

^ ^(2.-2,-,-,-)^... 



^ 5(V3/2) ^ 5(Vl/2) = ( 



12 2 2 
3' 9' 9' 9 



(49) 



Here the limiting Schmidt string Sq as in ()43|1 . is approached from above by 
iS(V„_i/2), and from below by iS(Vn_3/2), as n — oo. 

In this example, it seems that the notion of subspace entanglement has to 
be interpreted as a kind of mean entanglement per basis vector, rather than 
a total entanglement. Thus the 2-dimensional subspace V1/2, according to 
this notion, has a greater entanglement than, say, Vn~i/2i although the latter 
subspace may be of much greater dimension, containing many entangled 
states. 



6 Concluding remarks 

The notion of subspace entanglement that has been introduced has some 
interesting features. It is desirable in future work to try and establish that 
it does indeed make good sense in the context of applications to physics, 
perhaps along the lines that have been explored |ll]j in the case of state 
entanglement. To do that, it may well be necessary to relate more closely 
than we have done here, the entanglement of individual basis vectors in a 
subspace, with the notion of subspace entanglement. 

Are there other ways to measure subspace entanglement? Another way 
might be to consider an arbitrary orthonormal basis of the subspace, and to 
consider the system to be in a mixed state of those basis states, each with 
probability l/d, where d is the subspace dimension. Then we could associate 
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the entanglement of that mixed state with the entanglement of the subspace, 
using existing measures of entanglement of mixed states fTJ ^21 ^1 ■ The 
density operator for the mixed state in this case is simply a multiple (by 1/d) 
of the projection operator onto the subspace, so we would then be considering 
in a different way, the entanglement associated with a projection operator. 
It should be interesting to explore the connections between these two notions 
of subspace entanglement. 

We might also consider the possiblity that the subspace of states asso- 
ciated with a given quantum system is itself uncertain. In that situation 
it would seem appropriate to consider the extension of the superoperator 
density matrix M to the mixed case, with probabilties pi, p2, . . . , Pn associ- 
ated with different subspaces Vi, V2, . . . , Vn- Then we would need to extend 
existing notions of entanglement of mixed states to this new situation. 

We have considered above the coupling of an angular momentum j with 
a spin 1/2. There is a challenge to calculate the entanglement of irreducible 
subspaces with definite total angular momentum, in the case when two ar- 
bitrary angular momenta are coupled together. When this is done, it should 
be possible to see how the entanglement is related to the values of Clebsch- 
Gordan coefficients, in particular. This 'reduction entanglement' problem 
has an obvious extension to representations of other groups and algebras, 
and seems to open up a new aspect of the tensor product reduction problem 
in general, including cases involving infinite-dimensional representations, and 
cases involving tensor products of more than two representations. 

The concept of entanglement of a vector subspace of a tensor product 
space seems clearly to be of mathematical interest. It is less clear what 
may be its importance for physics, but we hope that the examples above are 
suggestive of important applications. 
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